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O I Abstract 

(D . 

Gauge fixed domain wall fermions are investigated in the reduced model at small Yukawa couplings. We present 
chiral propagators at the waveguide boundaries using quenched numerical simulations and analytic methods. 
O^ There is no evidence of mirror chiral modes at the waveguide boundaries. 
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^ ! Introduction 
> 

, I Recently it has been shown j^, || that nonpcrturbative gauge fixing Q can be apphed successfully to decouple the 
f— ^ . longitudinal gauge degrees of freedom from Abelian lattice chiral gauge theories in the limit of zero gauge coupling. 
[~^ ' After one gauge transforms a gauge-noninvariant lattice chiral gauge theory proposal like the Smit-Swift model or 
^D ■ the domain wall waveguide model, one picks up the longitudinal gauge degrees of freedom (radially frozen scalars) 
explicitly in the action. The job of gauge fixing is to find a phase transition where the gauge symmetry would be 
restored with the scalars decoupled. 

jrt I The gauge fixing approach deals with anomaly-free chiral gauge theories. It is manifestly local. The gauge fixing 

action involves a kinetic term for the longitudinal gauge degrees of freedom and allows for a renormalizable weak 

pi i' coupling expansion. 

Qj . The so called reduced model is defined by taking the limit of zero gauge coupling, after gauge transforming the 

gauge-noninvariant theory. Hence it has coupling between fermions and the scalars. The spectrum of the reduced 
model derived from the gauge fixed theory should be devoid of the scalars and consist only of free fermions in the 
appropriate representation of the chiral group. This is precisely what was achieved with recent investigations of 
Smit-Swift PI and domain wall models |2| and tuning counterterms posed no practical problems. In the reduced 

C^ ' model there is no constraint on the fermion representation from anomalies since the anomaly vanishes trivially in 
this case, and one can still study the fermion spectra. 

In this communication we extend our detailed study of the U(l) lattice chiral gauge theory with domain wall 
fermions and gauge fixing in the reduced model limit. The reduced model has a Yukawa coupling y which was put 
in by hand. Our previous study dealt with the situation at y — I and showed that the scalar fields were completely 
decoupled from the fermions and chiral modes were obtained only at the domain wall and at the anti-domain wall 
as in free domain wall fermions. This makes the gauge fixed domain wall model suitable for a U(l) chiral gauge 
theory at y — 1. On the other hand, at y = 0, fermion current considerations and also numerical simulations show 
that mirror chiral modes develop at the so-called waveguide boundaries making it unsuitable for the construction 
of a chiral gauge theory. Hence it is interesting to ask what happens at small nonzero values of Yukawa coupling. 
In general it is desirable to understand the properties of a model for all possible values of the coupling parameters. 
Here we investigate the gauge fixed domain wall waveguide model in the reduced limit for small nonzero values of 
the Yukawa coupling. 
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Gauge fixed Domain Wall Action 

Kaplan's free domain wall fermions [0 on a 4 + 1-dimensional L^Lg lattice (0 < s < L^ — 1, s labeling the 5-th 
dimension) with periodic boundary conditions in the 5-th direction and the domain wall mass taken as 

-mo, < s < Ls/2 
m{s)^ 0, s = 0,Ls/2 (1) 
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possess a lefthanded (LH) chiral mode bound to the domain wall at s = and a righthanded (RH) chiral mode 
bound to the anti-domain wall at s = Ls/2. For m^Ls ^ 1, these modes have exponentially small overlap. 

A 4-dimensional gauge field which is independent of s is then coupled to fermions only for a restricted number 
of s-slices around say, the anti-domain wall g, g] with a view to coupling only to the RH mode at the anti-domain 
wall. The gauge field is thus confined within a waveguide^ 

WG^{s:sq<s<si). (2) 

The waveguide boundaries sq and si should be reasonably far away from the domain wall and the anti-domain wall. 
Obviously, the hopping terms from sq to sq -I- 1 and that from si to si -I- 1 would break the local gauge invariance 
of the action. This is taken care of by gauge transforming the action and thereby picking up the longitudinal gauge 
degrees of freedom or radially frozen scalar fields 99 at the waveguide boundary, leading to the action (lattice constant 
is taken to be unity throughout): 
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where we have suppressed all indices other than s. tp and ip are the fermion fields, the projector Pl(r) is (1 T 75)/2 
and y is the Yukawa coupling introduced by hand at the waveguide boundaries. l/> [U) and W{U) are respectively 
the gauge covariant Dirac operator and the Wilson term (with Wilson r = 1) in 4 space-time dimensions. ^ and w 
are the 4-dimensional free Dirac and Wilson operators. 

The gauge fixed pure gauge action for C/(l), where the ghosts are free and decoupled, is: 

Sb{U) - Sg{U) + SgjiU) + S,t{U) (4) 

where, Sg is the usual Wilson plaquette action and Set is a counterterm. Previous studies 0, 13, R] show that only a 
gauge field mass counterterm is needed. Sgf is the gauge fixing term which is not just a naive lattice transcription of 
the continuum covariant gauge fixing, it has in addition appropriate irrelevant terms. As a result, Sgf has a unique 
absolute minimum at Uf^x — 1 validating weak coupling perturbation theory around gauge coupling g — Q. For an 
explicit expression of Sgf , see |3| . 

Obviously, the action Sb{U) is not gauge invariant. By giving it a gauge transformation the resulting action is 
S B {'PxUfj,x'Px+fi) ■ By restricting to the trivial orbit, U/^ix = 1, we arrive at the reduced model action 

Sreduced = Sf{U = 1) + SBivl 1 'Px+fi) (5) 

where Sp{U = 1) is obtained quite easily from eq.(w) and 

X X 

now is a higher-derivative scalar field theory action and Bx is given by, 

i?,^Vf ^^-^ + ^^- y with V,x=lTnrx^x+f.. (7) 



Perturbation theory around g — Q translates in the reduced model to the same around k — oo. 

In the following, we investigate the action (pi) at small nonzero y by numerical methods. To follow up the numerical 
investigation, analytic studies will then be made again at small y with slightly different boundary conditions. 
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Figure 1 : RR propagator at waveguide boundary s = 5 and LL propagator at waveguide boundary s 
and R = 0.2 {Ls = 22; a.p.b.c. in L4), (a) y = 0.25, (b) y = 0.50, (c) y = 0.75, (d) y = 1.0. 
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Numerical results at small Yukawa coupling 

At 2/ = 0, the domain wall and the anti-domain wall are detached from each other. In fact, in this case there is no 
coupling between the 4-dimensional worlds living on the So-th and the (sq + l)-th slices (and similarly for the Si-th 
and the (si + l)-th slices). From fermion current considerations, the LH chiral mode bound to the domain wall at 
s — Q would then necessitate the generation of a RH mode at s — sq and also at s = si + 1. Similarly a LH mode, 
mirror to the RH chiral mode at the anti-domain wall at s = Ls/'^^, would be generated at s = sq -I- 1 and also at 
s = si. Our numerical simulations at y = and similar studies in 0] clearly show that mirror chiral modes form at 
the waveguide boundaries. Generation of these mirror modes at the waveguide boundaries is independent of gauge 
fixing. 

On the other hand, with gauge fixing aX, y = 1 yj , the mirror chiral modes are certainly absent at the waveguide 
boundaries. In this case, the only chiral modes are at the domain wall and at the anti-domain wall and the spectrum 
is that of a free domain wall fermion. 

To investigate the interesting question what happens at small y, we looked for chiral modes at y = 0.75, 0.5, 0.25, 
especially at the waveguide boundaries. Specifically we numerically evaluate the LL and the RR propagators in 
momentum space conjugate to the 4 space-time dimensions and in configuration space for the 5-th (flavor-like) 
dimension. This is done using quenched configurations of the scalar fields at k = 0.05 and k — 0.2. This choice of 
the coupling parameters puts the theory quite close to the phase transition from a ferromagnetic broken phase (FM) 
to a rotationally noninvariant broken phase (FMD), staying within the FM phase (in the full theory with dynamical 
gauge fields FM-FMD is the transition where gauge symmetry is restored). For a detailed discussion of the phase 
diagram of the quenched reduced model, please see H. The (k, h) point we have chosen for the present investigation 
is also where we had previously done all our work at y = 1. 

On a 6^ X 16 lattice with L^ = 22 and ttiq = 0.5, the chiral propagators at the waveguide boundaries showed an 
increasing trend as the fourth component of momentum p4 (with p — Q) was decreased to the minimum value possible 
for this lattice size, something that could signify a pole at zero momentum, thereby indicating possible presence of 
mirror modes. The lower the Yukawa coupling, the sharper was the increasing trend of the chiral propagators with 
decreasing p^ . 

For a closer scrutiny we then took lattice sizes which were bigger in the 4-th direction to accommodate lower 
P4-values. We found that for each y there is a big enough lattice size for which the chiral propagators ultimately start 
showing a decreasing trend, indicating possible absence of mirror modes. Obviously the smaller the Yukawa coupling, 
the bigger the L4 extension was required. Moreover, all the chiral propagators at these small Yukawa couplings 
matched exactly with the corresponding free chiral propagators. These free chiral propagators were numerically 
determined by inverting the free fermionic action (obtained by putting 1^9 = 1 in the reduced fermionic action 
Sf{U = 1)) with the respective values of the Yukawa couplings (the free fermionic action at y = 1 is the free domain 
wall action of Kaplan Q). 



Our results at the waveguide boundaries s = sq = 5, and s = sq + 1 = 6 are summarized in Fig.n^ for the RR 
and the LL propagators respectively. In the figures, SF and FF respectively indicate data obtained by numerical 
simulation of the scalar- ferniion reduced model, eq. (o) , and by direct numerical inversion of the free fermion matrix 
(i.e. with if = 1) a,t the given Yukawa couplings. Error bars are smaller than the symbols. Dotted lines joining the 
data points are to guide the eye. Figure ^ also contains our previous data at y = 1 [0 for comparison. From the 
figures, existence of poles at zero 4-momentum for these chiral propagators does not seem likely. 

Similar numerical investigation is also carried out at the other waveguide boundary si, Si + 1. There too evidence 
for mirror modes is dim. 

Perturbation Theory at y ^ 1 

To confirm the indication from our numerical simulation that, for small y, the chiral propagators at the waveguide 
boundaries do not seem to have poles at zero momentum, it would be nice to calculate the propagators analytically. 
This can be done in perturbation theory in the coupling l/y/2K. This was also done for the fermion propagators to 
1-loop in our previous investigation at y = 1 . 

The present case, however, is more complicated due to the presence of two more defects created at the two 
waveguide boundaries because oi y ^ 1, in addition to the existing defects, namely, the domain wall at s = and 
the anti-domain wall at s = Ls/2. The defects are the places across which translational invariance is not maintained 
and boundary conditions need to be imposed to match the propagators from two sides of each defect. 

In this section we leave Kaplan boundary conditions and use Shamir boundary conditions S because imposing 
boundary conditions in the Shamir case would be a lot less tedious and our conclusions would qualitatively be the 
same as in the Kaplan case. 

The free Dirac operator in the Shamir case is defined on 4 space-time dimensions plus a finite 5-th dimension, 
< s < Lg — 1 with rigid walls at s = and s = Ls — 1. The domain wall mass is taken as (— ttiq) for all s— slices. 
Fermionic spectrum then consists of a LH chiral mode at s = and a RH chiral mode at s = Ls/2. (For further 
discussion see g and |^.) 

The waveguide in this case can be implemented by putting the same 4-dimensional gauge field over s-slices from 
So + 1 to Ls — 1 with a view to gauging only the RH mode at s = Ls — 1. The hopping terms from sq to sq -I- 1 
again break the local gauge invariance of the action. Longitudinal gauge degrees of freedom or radially frozen scalar 
fields if enter the action explicitly after a gauge transformation. The reduced model is then obtained by imposing 
U^. = 1: 



- y (^^""if^PLr''+'+T'^'ipPB.r'') ■ (8) 



The fermion propagators are obtained in momentum space for 4 space-time dimensions while staying in the 
coordinate space for the 5-th dimension. The procedure in the following is applicable to all values of the Yukawa 
coupling y including y = and y = 1. Needless to say, it is also applicable to Kaplan boundary conditions. 

In order to develop perturbation theory, in reduced model, we expand, 

^, = exp{ibe,) = i + ibe,-h^el + oib''), b = ^. (9) 

Thus the fermion action at the tree level is written as, 



^- ^ = E ^p ['^^-'t + ^^'^-y^st Pl + {Mpj ^^ Pr] V-* (10) 
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where p = sin(pp) and /5 = 7pP^, and 



(m|)^^ = Ml+y5s^s,+i5t,so (11) 



where y = 1 — y. The M and Af ^ are defined as, 

Mst = [l-mo + ^(l-cosp,0],5. 



M. 



.ft _ 






(12) 



The free fermion propagator can now formally be written as, 
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General solution of Gl is obtained by writing (|lj) explicitly: 

and similarly for Gr. We show only the calculations for obtaining Gl and henceforth drop the subscript L. 
We now define two translationally invariant regions / and // and set appropriate notation for G, 
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< s < So, 

So + I < s < Ls - 1, 
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G'^^^ and G*^^^ are now forced to satisfy the translationally invariant equation (with 5 = 1 — mp + X]u(l ~ cosp^)), 
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The above equation dl3) actually takes the system beyond the respective regions / and //. This will be taken care 
of by boundary conditions later. 

The solutions of cq.(n9) are expressed as sum of homogeneous and inhomogeneous solutions. 



G^^\p)s.t = g^^\t)e-'^' + h^^^\t)e"' + 
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where I = Ls/2, and 



cosh(a) = i[B + i^ 
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The third terms in eqs.(EG, |2l| ) are the inhomogeneous solutions, which are the same in both the regions. Avoiding 
singularity in a when B is zero restricts the allowed range of rriQ to < mp < 1. In this work we have taken toq = 0.5. 
The unknown functions g^^''^'{t) and h^^''^'{t) are determined from the following boundary conditions: 
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Figure 2: Tree level RR propagator at waveguide boundary s = 49 (512'^ x 8192 lattice, Lg — 100; a.p.b.c. in L4) 



where T — T — y{2 — y) and !F = p^ + I + B^ . Substituting 2i?sinli(Q:) sinli(a/) = X ^ and using the boundary 
conditions, eqs.(p3), we arrive at, 

A.v(i)=X(i), (24) 

where, v(t) = (5^-'^' ft,^^^ g^^^ ft,*^^^) is a 4-component vector, X(i) is another 4-coinponent vector and A is a 4 x 4 
matrix as given below, 
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The solution to the eqs.(p4|) is very complicated in general, particularly for finite Ls- However, g^^'^-'(i) and 
h^^''^\t) can be obtained for finite Lg by solving the above equations numerically for different t values. This way we 
can easily construct the tree level fermion propagators at any given s-slice. 

The solutions for (G'ij)^,* and the resulting propagators are obtained in exactly the same way. 

1-loop corrections to the tree level LL and RR propagators above can be easily found following and as in 
y = I are found to be negligible. This is because the scalar fields are almost completely decoupled. 

Using the analytic method described above, on a 512'^ x 8192 lattice with Lg = 100 and ttiq = 0.5, we plot in 
Fig.||tree level RR propagators for y = 0.25, 0.50, 0.75, 1.0 at the waveguide boundary s = sq — 49. From the figure 
it is clearly seen that there are no poles in these propagators, ruling out any possibility of a mirror mode at the 
waveguide boundary at any nonzero y (including arbitrarily small y). 



Discussion 

In our previous study |S| we investigated the U(l) gauge fixed domain wall waveguide model in the reduced limit 
at Yukawa coupling y — 1. Although the theory contained scalars explicitly in the action as a result of gauge 
noninvariance, we found no evidence of them. The fermionic spectrum was that of a free domain wall fermion. 

Since the Yukawa interaction couples only the two slices across a waveguide boundary, there is a new defect 
introduced at a waveguide boundary when y ^ 1. Especially at y — mirror chiral modes are developed at the 
waveguide boundaries. It is an interesting question to find out if these mirror modes persist to exist even for small 
values of the Yukawa coupling. 

Our numerical investigations (at small y) of chiral propagators at the waveguide boundaries initially show a 
disturbing trend of increasing with decreasing momentum. On bigger lattices and smaller momenta, however, this 
trend does not continue and poles of the chiral propagators do not seem likely. In addition, the data agree quite well 
with free chiral propagators obtained by numerically inverting the fermion matrix after putting ip = 1 for all the y 
values. Agreement with free propagators means that the scalar fields are decoupled. 

To be able to go to much smaller momentum, we then perform an analytic calculation of the propagators at y 7^ 1 
in the perturbative limit k — 00. Although done with different boundary conditions than we started with, qualitative 
conclusions can at least be made and the analytic tree propagators certainly show a very good resemblance to the 
numerically calculated ones and confirm the view that there are no poles for these propagators at the waveguide 
boundaries. 

Our conclusion would be that even for arbitrarily small or any nonzero Yukawa coupling, the fermion spectrum 
would be devoid of any ill-effects of the scalar fields. In fact the scalar fields are completely decoupled. This happens 
because the continuum limit is taken at the FM-FMD transition where the scalars do not scale, being far away 
from a ferromagnetic-paramagnetic transition. Obviously from a quenched numerical simulation alone, we cannot 
come to a firm conclusion. However, the 1-loop corrections to the perturbative tree level fermion propagators are 
negligibly small and fermion loops enter the calculation at least at the two loop level. Hence our conclusions about 
nonexistence of the mirror modes at the waveguide boundaries for small y seem to be correct even for the situation 
with dynamical fermions as long as the phase diagram stays qualitatively the same with dynamical fermions. Once 
the continuum limit can be taken at the FM-FMD transition, even an arbitrrily small y (which effectively mimics a 
scaling radial mode of the scalar fields) cannot produce the mirror modes and spoil the theory. 
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